Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 17, Issue 2, ICGST, Delaware, USA, December 2017 



ica«T 


www.icgst.com 


AC'-f'E 


Advanced D-Partitioning Stability Analysis of Digital Control Systems with 

Multivariable Parameters 


Kamen M. Yanev 1 

Department of Electrical Engineering, University of Botswana, Gaborone, Botswana 
yanevkm@yahoo.com, yanevkm@mopipi.ub.bw 
http ://www.ub .bw/ 


Abstract 


1. Introduction 


This paper contributes in further advancement of the 
method of D-Partitioning stability analysis applied to 
digital control system with multivariable parameters. It 
also explores the effects of simultaneous system 
uncertainties by determining graphically regions of 
stability in the space of the system’s parameters. The 
Euler's approximation is taken into account and the 
Bilinear Tustin Transform method is employed as one of 
the most precise facilities for systems discretization. The 
interaction between the variable parameters will bring a 
new light in the graphical solution of the problem of 
control system stability, where the analysis of control 
systems with two and three simultaneously variable 
parameters is taken into consideration. The proposed 
research is a continuation and further development of the 
already published author’s work on the advancement of 
the method of the D-Partitioning stability analysis of 
digital control systems with variable parameters. The 
suggested implementation for system’s stability analysis 
is fundamental for the further advancement of control 
theory in this area. 

Keywords: Multivariable Parameters, Stability Analysis 
Tool, D-Partitioning, Interaction, Regions of Stability; 


Nomenclature 

Go(s) Open-loop transfer function of the system 

G(s) Characteristic equation of the continuous system 

s Laplace operator 

K Variable gain of the control system 

Ti t T 2 T 3 Variable time constants of the control system 

T s Sampling period 

Tmin Minimum time-constant of the system 

cos Sampling frequency 

z z-transform operator 

ju, y System’ s variable parameters 

2- D Two-dimensional 

3- D Three-dimensional 

GM Gain margin 

PM Phase margin 


The D-partitioning method enables a quick and 
convenient determination of the regions of stability in 
case of variation of system’s parameters. Following some 
initial ideas of Neimark [1], the method was better 
clarified and further developed by the author in previous 
published work [2], [3], [4], [5], [6], [7], now calcified as 
Advanced D-Partitioning. After its further expansion, the 
method became a powerful tool for system stability 
analysis, having a number of advantages compared to 
other well-known stability analysis methods. It has the 
benefit of a clear graphical display of the variation of 
each parameter and its effect on the system’s stability. 

The current research is extending the strategy for stability 
analysis of digital control systems with multivariable 
parameters, suggesting a technique for solving the 
problem by applying the Advanced D-Partitioning 
method in the discrete time-domain. The Advanced D- 
Partitioning analysis of digital control systems with 
multivariable parameters is achieved by introducing 
discretization based on the Bilinear Tustin Transform [8], 
while the systems’ sampling time is chosen, considering 
the Euler's approximation [9], [10]. 

Exploring control systems with simultaneous variation of 
multiple parameters was inspired by cases of abnormal 
behaving systems, being stable at lower parameter 
values, further becoming unstable at a specific parameter 
range and turn out to be again stable again at higher 
parameter values [11], [12]. 

For the clarity of the presented research, this paper is 
organized in the following sequence. Initially, the nature 
of the D-Partitioning as stated by Neimark is explained. 
This is followed by further upgrading of the Advanced 
D-Partitioning analysis suggested by the author in his 
previous published work, now extended in the discrete 
time-domain. The results, achieved from the Advanced 
D-Partitioning analysis in the discrete-time domain for 
the cases of one, two and three variable parameters, are 
confirmed by the Nyquist and Bode stability criteria. The 
employed stability analysis methods are finally compared 
for the considered control systems variable parameters. 


9 



Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 17, Issue 2, ICGST, Delaware, USA, December 2017 


2. The Nature of the D-Partitioning as 
Originally Suggested by Neimark 

The position of the characteristic equation roots in the 
s-plane depends on the values of the system’s parameters. 
On this basis, Neimark suggested that the space of an 
n-order system’s characteristic equation coefficients can 
be partitioned into a number of regions corresponding to 
the number of roots in the left-hand side of the s-plane. 
This method was categorized as D-Partitioning. At its 
initial steps the method was stating rather its theoretical 
possibilities and its applications were quite limited. 
Although, Neimark expanded his initial ideas, they were 
rarely implemented because of their obscurity. 

3. Advanced D-Partitioning Stability Analysis 
(Case of One Variable Parameter) 

As already indicated, in previous published work, the 
author of this research suggested further advancement of 
the method of the D-Partitioning stability analysis of 
continuous control systems by one variable parameter. To 
facilitate and realize the stability analysis, the system’s 
characteristic equation is presented in the format [2], [3]: 

G(s) = a Q s n + a l s n 1 + ... + a n = 0 (1) 

Equation (1) can be also introduced in the following 
format to expose the variable parameter: 


The results obtained by the “dpartition” m-code are the 
same as those achieved by the “nyquist” m-code. This 
D-Partitioning procedure can work only on a computer 
where the new developed “dpartition” m-code is included 
in the MATLAB program. 

In order to benefit from the Advanced D-Partitioning 
analysis shown in this research, the wider engineering 
community can still use the “nyquist” m-code for the 
purpose of plotting the D-partitioning curve. 

The Advanced D-Partitioning analysis in case of one 
variable parameter is illustrated in this research by a 
system consisting of an armature-controlled dc motor and 
a type-driving mechanism is suggested [11], [12], [13]. 
Initially, it is considered that the system is experiencing 
variation of only one of its parameters. Usually, the 
variation of the load causes variation of the mechanism 
time-constant Ti. The transfer function of the open loop 
system in the continuous-time domain is presented as: 

10 

G 0 i ( s ) ~ — 

(1 + T lS )(\ + 0.5.s)(l + 0.8.S) 


T (0.4s 3 + 1.3s 2 + s) + 0.4s 2 + 1.3s + 1 

The characteristic equation of the continuous stand-alone 
unity feedback system is: 


G(s) = P(s) + vQ(s) = 0 (2) 

The D-partitioning regions could be obtained by 
substituting s =jco. 

G(jco) = P(jco) + vQ(jco) = 0 (3) 

Therefore, the variable parameter can be presented as a 
complex number: 


v 


P(jco) 

Q(jco) 


X(a>) + jY(a>) 


(4) 


The D-partitioning regions are obtained graphically in the 
v-plane by allocating values of the frequency within the 
range -oo < to < +oo. 

Since, in the s-plane, the region of stability remains on 
the left-hand side of the plane, in the complex plane 
v = X(a>) + jY(co) , a region of stability remains always 
on the left-hand side of the D-Partitioning curve for a 
change of frequency from -oo to +qo. 

The D-partitioning curve in terms of one variable 
parameter can be plotted in the complex plane within the 
frequency range -oo < oo < +oo, facilitated by MATLAB 
the “nyquist” m-code. The procedure can work on any 
computer where the MATLAB program is installed. 

To avoid any misunderstanding or misinterpretation of 
the D-Partitioning procedure, the “nyquist” m-code is 
modified into a “dpartition” m-code with the aid of the 
MATLAB Editor and a proper formatting. The 
“dpartition” m-code will plot the curve of a specific 
system parameter in terms of the frequency variation 
from -oo to +oo. 


G(s) = T (0.4s 3 + 1 ,3s 2 + s) + 0.4s 2 + 1 ,3s + 1 1 = 0 (6) 

The regions of stability reflecting the variation of the 
system’s parameter Ti are determined by implementing 
the method of the Advanced D-Partitioning. The time- 
constant Ti is presented from equation (6) as follows: 

TT s 2 +(T +T )s + K + \ 

T(s) = — — 1 = 

' TT/ + (T i+ T 2 )s 2 + s ( 6 ) 

_ 0.4s 2 + 1.3s + 11 

0.4s 3 + 1.3s 2 +s ' 


The variable Tj(s) is initially introduced as a continuous- 
time function and next converted into its digital 
equivalent Ti(z), facilitated by the Tustin Transform [14], 
[15], [16]. In accordance with the Euler's approximation, 
observing that T s < (0.1 T min to 0.2 T m i n ), the sampling 
period is chosen as T s = 0.05sec, since the continuous 
system minimum time-constant is T min = 0.5 sec. Then, 
the D-Partitioning in terms of the variable time-constant 
T](z) is achieved in the discrete-time domain by the code: 


» TI = tf([-0.4 -1.3 -11], [0.4 1.3 1 0]) 
Transfer function: 

-0.4 s A 2- 1.3 s -11 


0.4 s A 3 + 1.3 s A 2 +s 
» Tld = c2d(Tl,0.05, 'tustin ) 

Transfer function: 

-0.02536 z A 3 + 0.02002 z A 2 + 0.02377 z- 0.02161 


z A 3 - 2.844 z A 2 + 2.694 z - 0.8499 
Sampling time: 0.05 
» dpartition(Tld) 
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From where: 

, -0.02536z 3 + 0.02002z 2 + 0.02377z- 0.02161 n\ 

T,(z) = 5 5 ' 

z 3 - 2.844s 2 + 2.694z - 0.8499 


D-Partitioning in Terms of the Time-Constant T1 



Figure 1. Advanced D-Partitioning Analysis of the Digital Control 
System in Terms of the Variable Time-Constant Ti 

As seen from Figure 1, the Advanced D-Partitioning 
determines four regions in the complex Tl -Plane: D1(0), 
D2(0), D(l) and D(2). Only D1(0) and D2(0) are regions 
of stability, since the position of these regions is always 
on the left-hand side of the D-Partitioning curve within 
the range of co = ± co/2 = ±ln/lT s , = ± 62.83 rad/sec 
[17]. Only the real values of 77 are considered hence the 
regions of stability are reduced to lines of stability. 

4. System Evaluation with Nyquist and Bode 
Stability Criteria (Case of One Variable) 

If the time-constant is within 0 sec < 77 < 0.264 sec, 
corresponding to the region of stability D1(0), the system 
will be stable. If 77 is varied within 0.264 sec < 77 < 1 .48 
sec, the system will be unstable, being within the region 
of instability D(2). If 77 > 1.48 sec, the system is again 
stable and is operating in the region of stability D2(0). If 
the time-constant is 77 = 0.264sec or 77 = 1.48sec, the 
digital closed-loop control system becomes marginal. 

The open-loop transfer functions for both marginal cases 
of the control system in the continuous time-domain are 
determined accordingly: 

o 0.264 v 01 (1 + 0.264s)(l + 0.5s)(l + 0.8s) ^ 

10 

0. 1 056s 3 + 0.7432s 2 + 1 .564s + 1 




10 


01.48 or (l + 1.48s)(l + 0.5s)(l + 0.8s) 
10 


( 9 ) 


0.592s 3 + 2.324s 2 + 2.78s + 1 


The Advanced D-Partitioning in the discrete-time domain 
is confirmed with the Nyquist stability criterion and the 
Bode stability criterion for the two marginal values of 77. 
They are introduced into the system’s continuous transfer 
functions that are further converted into their digital 
equivalents [8], [13], [14], [17] with the Bilinear Tustin 
Transform. This is facilitated by the following code: 


» Gol=tf([0 10], [0.1056 0. 7432 1.564 1]) 
Transfer function: 

10 


0.1056 s A 3 + 0. 7432 s A 2 + 1.564 s + 1 
» Godl = c2d(Gol,0.05, ’ tustin ) 

Transfer function: 

0.001248 z A 3 + 0.003745 z A 2 + 0.003745 z + 0.001248 


z A 3 - 2.671 z A 2 + 2.375 z - 0. 7029 
Sampling time: 0.05 
» Go2=tf([0 10], [0.592 2.324 2. 78 1]) 
Transfer function: 

10 


0.592 s A 3 + 2.324 s A 2 + 2.78 s + 1 
» God2 = c2d(Go2,0.05, ’tustin') 

Transfer function: 

0.0002397 z A 3 + 0.0007191 z A 2 + 0.0007191 z + 0.0002397 


z A 3 - 2.811 z A 2 + 2.633 z - 0.8217 
Sampling time: 0.05 
» nyquist(Godl,God2) 

» bode(Godl,God2) 

» margin(Gpodl) 

GM = 0.00652 dB, PM = 0.0193 deg 
» margin(Gpod2) 

GM = - 0.0756 dB, PM = - 0.276 deg 


As seen, after applying the Bilinear Tustin Transform, the 
digital equivalents of equations (8) and (9) are reflected 
in equations (10) and (1 1), as seen from the code above. 

G (z) = G (z) = 

Od 0.264 od 1 

0.001248Z 3 + 0.003745Z 2 + 0,003745z + 0,001248 ^ 

z 3 -2.671z 2 +2.375z- 0.7029 
G (z) = G (z) = 

Od 1.48 od2 

3 2 (11) 

0.0002397Z 3 +0.0007191Z 2 +0.0007191Z + 0.0002397 v 2 
z 3 -2.811z 2 +2.633z-0.8217 

As a result of the applied code, Figure 2 demonstrates the 
Nyquist stability criterion in the discrete-time domain for 
both of the cases related to marginal time-constant values 
77 = 0.264sec and 77 = 1 .48sec. 


Marginal Nyquist Diagrams in the Discrete Time-Domain 
Cases of Tl = 0.264 sec and Tl = 1.48 sec 


Marginal System 
. Tl = 0.264 sec 

T1= 1.48 sec / 

+62.83 
lad sec 

System: God2 

X Real: -1 

f Imag: 0.00609 

1 Frequency (rad/sec): -2.26 


-62.83 

rad/sec 


■ 6 1 1 ‘ 1 — 

-4 -3 -2 -1 0 1 

Real Axis 


Figure 2. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Nyquist Stability Criterion in the 
Discrete-Time Domain (Marginal Cases Ti = 0.264sec; Ti = 1.48sec) 
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As seen from Figure 2, the Nyquist curves of the digital 
open loop system coincide and are passing approximately 
via the point (-1, jO), hence the closed loop digital control 
system is marginal. This confirms the results, obtained 
from the Advanced D-Partitioning stability analysis. 

The achieved results from the code, confirm also the 
outcome of the Advanced D-Partitioning analysis with 
the aid of the Bode stability criterion in the discrete-time 
domain for both of the cases related to marginal time- 
constant values Ti = 0.264sec and Ti= 1.48sec. 

Marginal Bode Diagrams in the Didcrete Time-Domain 
Cases of T1 = 0.264 sec and T1 = 1.48 sec 



Frequency (rad/sec) 

Figure 3. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Bode Stability Criterion in the 
Discrete-Time Domain (Marginal Cases Ti = 0.264sec; Ti = 1.48sec) 

As shown in Figure 3, if the marginal value of the time- 
constant is Tj = 0.264 sec, the results for the digital 
system gain margin is GM = 0.00652 dB « 0 dB and the 
system phase margin is PM = 0.0193° « 0°. The results 
also prove the second marginal value of the time-constant 
Ti = 1.48sec, where the gain margin is GM = - 0.0756 
dB « 0 dB and phase margin is PM = - 0.276° « 0°. 

Again, the analysis of both Nyquist stability criterion and 
Bode stability criterion is restricted to the frequency 
co — i co s /2 = + 2n/2T s = ± 62.83 rad/sec [11], [12], [17]. 

The system stability within the regions D1(0) and D2(0) 
is explored in the discrete-time domain by allocating 
values for the variable time-constant Ti = 0.1 sec and 
Ti = 5 sec and applying the following code: 


» Go3=tf([0 10], [0.04 0.53 1.4 1J) 

» God3 = c2d(Go3,0.05, ' tustin ) 

Transfer function: 

0.002886 z A 3 + 0.008658 z A 2 + 0.008658 z + 0.002886 

z A 3 - 2.444 z A 2 + 1.956 z - 0.51 
Sampling time: 0.05 
» Go4=tf([0 10], [2 6.9 6.3 1]) 

» God4 = c2d(Go4,0.05, 'tustin ) 

Transfer function: 

7.1 79e-005 z A 3 + 0.0002154 z A 2 + 0.0002154 z + 7.179e-005 

z A 3 - 2.834 z A 2 + 2.676 z - 0.8415 
Sampling time: 0.05 
» nyquist(God3,God4) 

» bode(God3,God4) 

» margin(God3) 

GM = 4.89 dB, PM = 15.5 deg 
» m argin (God4) 

GM = 6.34 dB, PM = 24 deg 


Nyquist Diagrams in the Discrete Time-Domain 
Cases of Stable System: TI = 0.1 sec and TI = 5 sec 



Figure 4. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Nyquist Stability Criterion in the 
Discrete-Time Domain (Stable Case Ti = O.lsec; Ti = 5sec) 

Since, the Nyquist curves of the digital open loop control 
system are passing approximately via the points at the 
real axes (-0.566, jO) and (-0.474, jO), the closed loop 
digital control system is stable for both cases, confirming 
the result from the Advanced D-Partitioning analysis. 


Bode Diagram in the Discrete Time-Domain (Case TI = 0.1 sec) 
Gm = 4.89 dB (at 5.87 rad/sec) , Pm = 15.5 deg (at 4.46 rad/sec) 



io 2 io 1 10 ° io 1 io 2 

Frequency (rad/sec) 


Figure 5. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Bode Stability Criterion in the 
Discrete-Time domain (Case of Stable System Ti = O.lsec) 


Bode Diagram in the Discrete-Time (Case TI = 5 sec) 

Gm = 6.34 dB (at 1 .77 rad/sec) , Pm = 24 deg (at 1 .21 rad/sec) 



Figure 6. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Bode Stability Criterion in the 
Discrete-Time domain (Case of Stable System Ti = 5sec) 
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As seen from Figure 5 and Figure 6, and from the results 
of the code, the Advanced D-Partitioning analysis is 
confirmed with the aid of the Bode stability criterion in 
the discrete-time domain for both of the cases related to a 
stable system ( Ti = 0.1 sec and Ti = 5sec). The positive 
gain and phase margins ( GM= 4.89 dB and PM = 15.5 °) 
and ( GM= 6.34 dB and PM= 24 °), prove that the digital 
system is stable within the regions D1(0) and D2(0). 

Further, the system stability within the region D(2) is 
explored in the discrete-time domain by allocating a 
value for the variable Ti = 1 sec and applying the code: 


» Go5=tf([0 10], [0.4 1.72.3 1]) 

» God5 = c2d(Go5,0.05, ’ tustin ) 

Transfer function: 

0.000352 z A 3 + 0.001056 z A 2 + 0.001056 z + 0.000352 


z A 3 - 2. 795 z A 2 + 2. 604 z - 0. 8085 
Sampling time: 0.05 
» nyquist(God5) 

» bode (God 5) 

» margin(God5) 

GM= - 1.14 dB, PM= - 4.11deg 


Nyquist Diagram in the Discrete time-Domain 
Case of Unstable System (TI = 1 sec) 



If applying the Nyquist stability criterion for the case 
Ti = 1, the Nyquist curve of the digital open loop control 
system is passing approximately via the points at the real 
axes (-1.14, jO). Therefore, the closed loop digital control 
system is unstable. 

By applying the Bode stability criterion for the case 
Ti = 1, the negative values of the digital system margins 
( GM = -1.14 dB and PM = - 4.11 °), confirm that the 
digital system is unstable within the region D(2). 

In case of control systems with one variable parameter, 
the Advanced D-Partitioning method directly exposes the 
transparent images of regions of stability and instability, 
as well as the system’s margins of stability in the 
complex plane of the variable parameter. It is also 
demonstrated that the Advanced D-Partitioning can be 
successfully applied to digital control systems. Both 
Nyquist and Bode stability criteria confirm the results 
from the Advanced D-Partitioning analysis. 

5. Advanced D-Partitioning Stability Analysis 
(Case of Two Variable Parameters) 

Further upgrade of the Advanced D-Partitioning analysis 
in case of systems with two simultaneously variable 
parameters [5], [6], [7] is suggested in this research. The 
system’s characteristic equation generally is presented as: 

G(s ) = juP(s) + yQ(s) + R(s) = 0 (12) 

where P(s), Q(s), and R(s) are polynomials of s 
ju and yare system’s variables parameters 

The border of the D-Partitioning regions in the plain 
(//, y) is determined by: 

G(j(D) = juP(jcD) + yQ(jco) + R(jco) = 0 (13) 

If a hypothetical third order unity feedback system of 
Type 0 is considered, the system’s characteristic equation 
can be presented as: 

G(s) = (T x s + 1 )(T 2 s + l)(r 3 s + 1) + K = 0 (14) 


Figure 7. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Nyquist Stability Criterion in the 
Discrete-Time domain (Unstable Case Ti = 1 sec) 


Bode Diagram in the Discrete Time-Domain (Case TI = 1 sec) 
Gm = -1 .14 dB (at 2.4 rad/sec) , Pm = -4.1 1 deg (at 2.54 rad/sec) 



It is suggested that two of the system’s parameters are 
simultaneously variable: 

T! = T=jU , K=y. (15) 

The objective is to determine the regions of variation of 
these two parameters, for which the system will be stable. 
Equations (15) are substituted in (14), from where: 

u[T T s 3 + (T + T )s 2 + s]+ 1 

2 3 2 s' I (16) 

+ r + T T s 2 + (T + r )s + 1 = 0 J 

2 3 2 3 


By substituting s = jco in equation (16) and taking into 
account equation (13), equation (16) could be presented 
in the detailed form as set of equations (17): 


P(jco) = [T T {jcof + (T + T + ja > ] 

2 3 2 3 

Q(ja) = 1 

R(jco) = TT(jco) 2 +(T + T )ja> + 1 

2 3 2 3 


(17) 


Figure 8. Confirmation of the Result of the Advanced D-Partitioning 
Analysis with the Aid of the Bode Stability Criterion in the 
Discrete-Time domain (Unstable Case Ti = 1 sec) 


Taking into consideration equations (13) and (17), the 
variable system’s parameters can be determined as: 
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r = 


A T +T 

_L _ 2 3 

A TTco 2 -\ 

2 3 

A (TT<o 2 -lf+(T +T) 2 co 2 +\ 
A 2 = 2 3 2 3 

A 


TTco - 1 

2 3 


or A = T 2 T 3 co 


(18) 


(19) 


The determinant A becomes A = 0 at a specific frequency 
co = (Ooo that can be found out from equation (19) as: 


1 

T T 

2 3 


( 20 ) 


If A = 0, it is clear from equations (20) that at co = C 0 oo 
both system parameters are approaching infinity: 

K® „)->°°, ( 2 i) 

This implies that the main D-Partitioning curve has an 
interruption, or a breakdown, at a frequency co = C0oo.lt 
consists of two parts, plotted within the frequency ranges 

0 < co < cOoo and co x < co < oo . 

For a better clarification and simplicity, first the 
functions p(co) and /(co) are plotted, as shown in Figure 
9(a) and Figure 9(b). 



L/ 


■ — kr 
1 \ “ 

1 


& = GO 


(b) 

Figure 9. The Graphical Presentations of jli(co) and y(ce>) Showing the 
Interruption of the Curves at a Frequency go = (Ooo 

The regions of the D-Partitioning also depend on straight 
lines in the (p ,y) plane, known as special lines. The 
special lines are plotted for two border frequencies co = 0 
and co = +oo. Then the coefficients a„ and a 0 of the 
equation ( 1 ) depend directly on the parameters p and y 
and the equations of the special lines are obtained by: 


a n = JL1T2T3 = 0 , 
jLi = 0 = Const 


a Q = y + 1 = 0 
y = - 1 = Const 


} 


( 22 ) 


By combining the curves ju(co), y(co) and the special lines, 
the D-Partitioning is obtained in the (p , y) plane as seen 
in Figure 10. Taking into account that p = 71 = T is a 
time-constant and it can adopt only positive values, only 
the stable region Di(0) should be considered. The region 
of stability Di(0) is locked within the left-hand side of the 
D-Partitioning curve, corresponding to frequency rise 
from co = Ox to co —> 00 and the special line y = -1. 
Since the gain K = y may also adopt only positive 
values, the realistic border of the stable region Di(0) 
should be considered K = y = 0. 



Figure 10. The Regions of Stability Di(0) and D 2 (0) Locked between 
the D-Partitioning Curves and the Special Lines 

If the gain K < y min , the system is stable for any values of 
the time-constant p = T\. For K > y min , the system is 
stable only for very small or very large values of the 
time-constant p = T\. This case is demonstrated for the 
system of the armature-controlled dc motor and a type- 
driving mechanism with open-loop transfer function: 

K 


G 02( S ) = - 


(23) 


(1 + Ts)(l + 0.5s)(l + 0.8s) 

The characteristic equation of the feedback system is: 

K + (1 + Ts)(l + 0.5s)(l + O. 85 ) = 0 (24) 

By substituting s =j equation (24) is modified to: 

K = -1 + (1.3T + 0 A)co 2 + (25) 

+ j co (OAT co 2 -l .3 - T) 

Since the gain may obtain only real values, the imaginary 
term of equation (25) is set to zero, from where: 

„ 1.3 + T 


OAT 


(26) 


The result of (26) is substituted into the real part of 
equation (25), from where: 


JC = 13P + 1.69T + a52 U 

0.4T T 


(27) 


When considering the continuous control system with the 
transfer function of equation (23), the D-Partitioning 
curve K = / (T) is plotted with the aid of the code: 


»T = 0:0. 1:5; 

»K = 3.25. *T+4.225+1.3./T 

K = 

Columns 1 through 10 

Inf 17.5500 11.3750 9.5333 8.7750 

8.4500 

8.3417 

8.3571 8.4500 8.5944 

Columns 11 through 20 

8.7750 8.9818 9.2083 9.4500 9.7036 

9.9667 

10.2375 

10.5147 10.7972 11.0842 

Columns 21 through 30 

11.3750 11.6690 11.9659 12.2652 

12.5667 

12.8700 

13.1750 13.4815 13.7893 14.0983 

Columns 31 through 40 

14.4083 14.7194 15.0313 15.3439 

15.6574 

15.9714 

16.2861 16.6014 16.9171 17.2333 

Columns 41 through 50 

17.5500 17.8671 18.1845 18.5023 

18.8205 

19.1389 

19.4576 19.7766 20.0958 20.4153 

Column 51 

20.7350 
» plot(T,K) 
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D-Partitioning by Two Variable Parameters (System Type 0) 



As seen from Figure 1 1, if K < 8.3417, a limit determined 
with the aid of MATLAB interactive procedure, the 
system is stable for any value of the T. The system 
performance and stability depends on the interaction 
between the two simultaneously varying parameters. 

Further, a very important phenomenon is demonstrated in 
Figure 1 1 . It is seen that the D-Partitioning curve K =f(T) 
defines the border between the region of stability D(0) 
and instability D(l) for the case of simultaneous variation 
of the two system parameters. Each point of the D- 
Partitioning curve represents the simultaneous 
marginal values of the two variable parameters. This 
result is a unique advancement and innovation in the 
theory of control systems stability analysis. 

The system performance is investigated for any values of 
the variable gain K, like K = 10, 12, K = 14. If K is 
varied, this affects the values of T at which the system 
may become unstable. Higher values of the gain K , 
enlarges the range of the time-constant T at which the 
system will fall into instability. 

The parameters interaction of the continuous control 
system is demonstrated in case of variation of the time- 
constant T, when the gain is set to K = 10. If the time- 
constant range is 0 < T < 0.25 sec, or T > 1.5 sec the 
system is stable. The system becomes unstable in the 
time-constant range 0.25 sec <T< 1.5 sec. 

If the results of the continuous control system are 
compared with the results for the same system in the 
digital time-domain, as demonstrated in Figure 1, it is 
obvious that there is a very close match. In the digital 
time-domain within the range of co= ± coj. 2 = ±2rc/2T s , 
if the gain is set to K = 10, the system is stable for the 
time intervals 0 < T < 0.264 sec, or T > 1.48 sec. It is 
unstable in the range 0.264 sec < T < 1.48 sec. This 
insignificant difference of the results is due to the 
implementation of the Euler's approximation. 

6. System Evaluation with Nyquist and Bode 
Stability Criteria (Case of Two Variables) 

The performance of the control system with two variable 
parameters can be tested with the Nyquist stability 
criterion and the Bode stability criterion in the discrete- 
time domain for any of the simultaneous marginal values 


of the variables T and K. Initially, the system’s transfer 
function is introduced as a continuous-time function and 
next it is converted into its digital equivalent. 

By substituting any two simultaneous marginal values of 
the two variable parameters, for instance, time-constant 
T= 0.25 sec and gain K = 10 in equation (23) results in: 


(28) 


Again, taking into account the Euler's approximation, 
T s < (0.1 T min to 0.2 Tmin), where the system’s minimum 
time-constant is T min = 0.5 sec and the sampling period is 
chosen as T s = 0.05sec. The system’s stability assessment 
in the discrete-time domain is achieved by the code: 


» Go21=tf([0 10], [0.1 0.725 1.55 1]) 

Transfer function: 

10 

0.1 s A 3 + 0. 725 s A 2 + 1.55 s + 1 
» God21 = c2d(Go21,0.05, 'tustin ) 

Transfer function: 

0.001312 z A 3 + 0.003935 z A 2 + 0.003935 z + 0.001312 

z A 3 - 2.662 z A 2 + 2.359 z - 0.6954 
Sampling time: 0.05 
» nyquist(God21) 

» m argin ( God21 ) 

GM = 0.017 dB, PM = 0. 028 deg 


G 0 (s) = 


10 


(1 + 0.25s)(l + 0.5s)(l + 0.8s) 
10 

n n 3 +n 7?Sc 2 +i ssc + i 


When applying the Nyquist stability criterion, as seen 
from Figure 12, the Nyquist curve of the digital open 
loop system is passing approximately via the point (-1, 
jO), therefore the closed loop digital control system is 
marginal. This confirms the results, obtained from the 
Advanced D-Partitioning stability analysis for the 
marginal system with two variable parameters. 


Prove of the D-Partitioning Analysis for the Case of Two 
Simultaneously Marginal Parameter Values T = 0.25 sec: K = 10 



Figure 12. Confirmation of the result of the Advanced D-Partitioning 
Analysis with the Aid of the Nyquist Stability Criterion in the 
Discrete-Time domain (Marginal Case T = 0.25sec; K = 10) 


As seen from the Bode diagram presented at Figure 13, 
the system is explored in the discrete-time domain for 
T= 0.25 sec and K = 10. The achieved results for the gain 
and phase margins are GM = 0.204 dB « 0 dB and 
PM= 0.729° « 0°, proving that the system is marginal. 
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Prove of the D-Partitioning for the Case of Two Marginal Parameters 
Gm = 0.204 dB (at 3.92 rad/sec) , Pm = 0.729 deg (at 3.88 rad/sec) 



Figure 13. Prove of the D-Partitioning Analysis with the aid of the 
Bode Stability Criterion - Case of Two Simultaneously Marginal 
Values (Case T = 0.25sec and K = 10) 

The system stability within the region D(0) is explored in 
the discrete-time domain if T = 1 sec and K = 5. The 
positive values for the gain and phase margins ( GM = 
4.89 dB and PM = 20.6°), prove that the digital system is 
stable within the region D(0), as seen from the code: 


» Go22=tf([0 5], [0.4 1 . 7 2.3 1]) 
Transfer function: 

5 


0.4 s A 3 + 1.7 s A 2 + 2.3 s + 1 
» God22 = c2d(Go22,0.05, Austin ) 

Transfer function: 

0.000176 z A 3 + 0.0005279 z A 2 + 0.0005279 z + 0.000176 


z A 3 - 2. 795 z A 2 + 2.604 z - 0.8085 
Sampling time: 0.05 
» margin(God22) 

GM= 4.89 dB, PM = 20.6 deg 


The system instability within the region D(l) is explored 
in the discrete -time if T= 1 sec and K = 14. The negative 
values of the digital system margins ( GM= -4.06 dB and 
PM = - 13.7°), prove that the digital system is unstable 
within the region D(l), as seen from the code: 


» Go23=tf([0 14], [0.4 1.72.3 1]) 
Transfer function: 

14 


0.4 s A 3 + 1.7 s A 2 + 2.3 s + 1 
» God23 = c2d(Go23,0.05, 'tustin ) 

Transfer function: 

0.0004927 z A 3 + 0.001478 z A 2 + 0.001478 z + 0.0004927 


z A 3 - 2. 795 z A 2 + 2.604 z - 0.8085 
Sampling time: 0.05 
» margin(God23) 

GM = -4.06 dB, PM = -13. 7 deg 


7. Advanced D-Partitioning Stability Analysis 
(Case of Three Variable Parameters) 


Further advancement of the D-Partitioning is suggested 
in case of systems with three simultaneously variable 
parameters, where their interaction is presented in the 
3 -Dimensional (3-D) space. This is demonstrated for the 
control system of the armature- controlled dc motor and a 
type-driving mechanism, where the system’s gain K and 



two of the time-constants Ti and T 2 are variable. The 
open-loop transfer function of the system is presented as: 

G(V) = — ( 29 ) 

( 1 + T iS )(\ + T 2 s)(1 + 0.&s ) 

The characteristic equation of the feedback system is: 

K + (1 + 7»(1 + T 2 s)(l + 0.8 s) = 0 (30) 

By substituting s =jco equation (30) is modified to: 

J A = -i + o.8(r 1 +r 2 )6f+ ( 31 ) 

+ j co{— T x — T 2 — 0.8 + 0. $T X T 2 ooP) 


The imaginary term of equation (31) is set to zero, since 
the gain may obtain only real values, from where: 


2 T x + T 2 + 0.8 
“ 0.ST x T 2 


(32) 


When the result of (32) is substituted into the real part of 
equation (31), results in: 


^ =-1+ 0.8(7; +r 2 )+T±Zi±M 
o.87;r 2 


(34) 


The D-Partitioning plane is created in the (K, Tj, Ti )- Space 
of the variable parameters. The code is only partly 
presented, due to its considerable length: 


» [X,Y] = meshgrid([0:.15:3]) 

X = 

Columns 1 through 10 

0 0.1500 0.3000 0.4500 0.6000 0.7500 0.9000 1.0500 

1.2000 1.3500 

0 0.1500 0.3000 0.4500 0.6000 0.7500 0.9000 1.0500 

1.2000 1.3500 


Columns 11 through 20 
1.5000 1.6500 1.8000 

1.9500 

2.1000 

2.2500 

2.4000 

2.5500 2.7000 2.8500 
1.5000 1.6500 1.8000 

1.9500 

2.1000 

2.2500 

2.4000 

2.5500 2.7000 2.8500 

Y = 

Columns 1 through 10 

0 0 0 0 

0 0 

0 

0 0 

0 

0.1500 0.1500 0.1500 

0.1500 

0.1500 

0.1500 

0.1500 

0.1500 0.1500 0.1500 
0.3000 0.3000 0.3000 

0.3000 

0.3000 

0.3000 

0.3000 

0.3000 0.3000 0.3000 






Column 21 
0 

0.1500 

0.3000 

»Z = -1+0.8. *(X+ Y) +(X+ Y+0. 8)./(0. 8. *X *Y) 
Warning: Divide by zero. 

Z = 

Columns 1 through 10 


Inf 

Inf 

Inf Inf Inf 

Inf 

Inf Inf 

Inf 

Inf 






Inf 

60.3511 

34.0822 25.4059 

21.1278 

18.6089 16.9696 

15.8330 

15.0106 

14.3975 




Inf 

34.0822 

18.9244 13.9519 

11.5256 

10.1178 

9.2193 

8.6117 

8.1861 

7.8817 





3.6249 3.5605 


Column 21 
Inf 

12.4922 

7.3344 


4.7444 

» surf(X, Y,Z,gradient(Z)) 
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D-Partitioning Analysis in Terms of the 
Three Simultaneously Variable Parameters 


* 

ii 

N 


Y = T2 



Figure 14. D-Partitioning in Case of a System with Three 
Simultaneously Variable Parameters 

The D-Partitioning plane, as seen in Figure 14, is 
defining the border between the space of stability D(0), 
below the plane, and the space of instability D(l), above 
the plane, for the case of simultaneous variation of three 
system’s parameters. 

According to the Advanced D-Partitioning rule, in this 
case the space of stability is the one that remains always 
on the left-hand side of the D-Partitioning plane within 
the range of co = ± co s f 2 = ±2n/2T s . Following this rule, 
the stable space is proved to be the space D(0). 

Any point of the D-partitioning plane represents the 
simultaneous marginal values of the three variable 
parameters, which is a unique advancement and in 
the theory of control systems stability analysis. 


8. System Evaluation with Nyquist and Bode 
Stability Criteria (Case of Three Variables) 


For the chosen point on the D-partitioning plane, shown 
in Figure 14, the system marginal gain is K = Z = 9.034 
and the marginal values of the two time-constants are 
Ti = X= 0.4141 sec and T 2 = 7= 0.6172 sec accordingly. 
To test the system’s performance, its transfer function is 
introduced initially in the continuous-time domain and 
then converted into its digital equivalent. Substituting the 
three marginal values of the variable parameters in 
equation (29) results in: 


G 03l( S ) = 


9.034 

(1 + 0.414 ls)(l + 0.61 72s)(l + 0. 8s) 
9.034 


0.192s 3 + 1.04s 2 +I.85 + I 



The system stability within the space D(0) is explored in 
the discrete-time domain if the following values are 
allocated to the variable parameters K = 3, Ti = 0.4 141 sec 
and T 2 = 0.6172sec, resulting in equation (34). 


^032 ( 5 ) — 


(1 + 0.4 14 ls)(l + 0.61 72s)(l + 0.8s) 
3 


0.192s 3 +1.04 s 2 +1.8s + 1 


(34) 


The system instability within the space D(l) is exploredin 
the discrete-time domain if K= 20, Tj = 0.4141sec and 
T 2 = 0.6172sec, resulting in equation (35). 


G Q33 (s) = 


20 

(1 + 0.4141s)(l + 0.6172s)(l + 0.8s) 
20 


0.192s 3 + 1. 04s 2 +1. 8s + 1 


(35) 


Then the system’s Nyquist and Bode stability assessment 
in the discrete-time domain is achieved by the code: 


» Go31=tf([0 9.034], [0.192 1.04 1.8 1]) 
Transfer function: 

9.034 


0.192 s A 3 + 1.04 s A 2 + 1.8 s + 1 
» God31 = c2d(Go31,0.05, Austin ) 

Transfer function: 

0.0006417 z A 3 + 0.001925 z A 2 + 0.001925 z + 0.0006417 


z A 3 - 2. 742 z A 2 + 2.505 z - 0. 7626 
» Go32=tf([0 3], [0.192 1.04 1.8 1 ]) 
Transfer function: 

3 


0.192 s A 3 + 1.04 s A 2 + 1.8 s + 1 
» God32 = c2d(Go32,0.05, Austin ) 

Transfer function: 

0.0002139 z A 3 + 0.0006417 z A 2 + 0.0006417 z + 0.0002139 


z A 3 - 2. 742 z A 2 + 2.505 z - 0. 7626 
» Go33=tf([0 20], [0.192 1.04 1.8 1]) 
Transfer function: 

20 


0.192 s A 3 + 1.04 s A 2 + 1.8 s + 1 
» God33 = c2d(Go33,0.05, 'tustin ) 

Transfer function: 

0.001426 z A 3 + 0.004278 z A 2 + 0.004278 z + 0.001426 


z A 3 - 2. 742 z A 2 + 2.505 z - 0. 7626 
Sampling time: 0.05 
» nyquist(God31,God32,God33) 

» margin(God31) 

GM = -0.224 dB, PM = -0.921 deg 
» margin(God32) 

GM = 9.3 dB, PM = 45 deg 
» margin(God33) 

GM = -7.18 dB, PM = -22.8 deg 


Nyquist Diagrams Proving the D-Partitioning for Cases of 
Marginal. Stable and Unstable System with Three Variables 



Figure 15. Prove of the D-Partitioning Analysis - Cases of Marginal, 
Stable and Unstable System with Three Variable Parameters 
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Applying the Nyquist stability criterion for the cases of 
marginal, stable and unstable system with three variable 
parameters, prove exactly the results from the Advanced 
D-Partitioning in the discrete time-domain. 

The system’s marginal case is also confirmed with the 
aid of the Bode stability criterion as shown in Figure 16. 

Bode Diagram Proving the Marginal Case 
Gm = -0.244 dB (at 3.06 rad/sec). Pm = -0.921 deg (at 3.1 rad/sec} 



Figure 16. Prove of the D-Partitioning Analysis - Case Marginal 
System (Case Ti = 0.4 141 sec, T 2 = 0.6172sec and K = 9.034) 

If the variable parameters are K = 9.034, Ti = 0.4 141 sec 
and T 2 = 0.6172sec, the results for the digital system gain 
margin is GM = —0.244 dB « 0 dB and the system phase 
margin is PM = —0.921° « 0°, as seen in the Bode 
diagrams in Figure 16. This proves that the three variable 
parameters have simultaneously marginal values. 

The system stability within the space D(0) is explored in 
the discrete-time domain if the following values are 
allocated to the variable parameters K = 3, Ti = 0.4 141 sec 
and T 2 = 0.6172sec. As seen in the code, the positive 
values of the gain and phase margins (GM =9.3 dB and 
PM = 45°) prove that if the system’s gain is reduced to 
K = 3, while keeping the values of the time constants, the 
digital control system will operate in the space of 
stability D(0). Similar result can be achieved if any of the 
two time constant values is reduced. 

The system instability within the space D(l) is explored 
in the discrete-time domain if K = 20, 77 = 0.4 141 sec 
and T 2 = 0.6172 sec. From the code above, the negative 
values of the gain and phase margins (GM = -7.18 dB 
and PM = -22.8°) prove that if the system’s gain is 
increased to K = 20, while keeping the values of the time 
constants, the digital control system operates in the space 
of instability D(l). Similar result can be achieved if any 
of the two time constant values is increased. 


9. Conclusion 


The main contribution of this research is the achieved 
further development of the Advanced D-Partitioning 
method in its application for the analysis of digital 
control systems with multivariable parameters. Cases of 
systems with one, two or three simultaneously variable 
parameters are illustrated. There is a close match between 
the discrete and the continuous-time system performance 
if the Euler's approximation is observed. 


The successful application of the analysis method in the 
discrete-time domain, by employing the Bilinear Tustin 


Transform, demonstrates that the basic concept of the 
Advanced D-Partitioning is relevant and valid for digital 
control systems as well. The control system’s 
discretization, based on the Bilinear Tustin Transform is 
providing the most precise results. 

When the Advanced D-Partitioning is compared with the 
Nyquist and the Bode stability criteria, even if both these 
methods are applicable for continuous and for digital 
control systems, their disadvantage is that they cannot 
directly establish the marginal value of a variable 
parameter. In this case, they can be used only to confirm 
the results from the Advanced D-Partitioning method. 

The research further contributes to knowledge, since the 
Advanced D-Partitioning is the only method of stability 
analysis that in terms of its results illustrates graphically 
the regions of stability and instability in case of system’s 
multivariable parameters. It was demonstrated that any 
point of the D-partitioning curve or D-partitioning plane, 
represents the simultaneous marginal values of the two or 
three variable parameters. This is a unique innovation in 
the theory of control systems stability analysis. 
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